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^ ■ Abstract. Correlation between emitted particles from a thermalized medium 

£h . carry important information about the space-time extent and the dynamics of 

the particle emitting source. Pion emission source functions Sir) have been 
measured in PHENIX using an imaging technique, and it has been found that 

■ it contains a heavy (power-law) tail. It leads us to the assumption that there 
0^ \ is a halo of secondary pions, surrounding the core system that undergoes the 

■ hydrodynamical evolution. THERMINATOR, a Monte Carlo event generator 
designed for studies in relativistic heavy-ion collisions is used to model and 
investigate the influence of resonance decays on the tail of particle emission 
source. It lacks the implementation of rescattering, which, as we have learnt 
from previous Hadron Rescattering Code simulations, is a possible explanation 
of the power-law tail. Our studies also have shown that none of the pion sources 
alone are responsible for the tail: either the resonance decays or rescattering 
or a combination of these effects is a possible cause of the heavy tail. 
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1. Introduction 

Correlation between identical pions emitted from a thermalized source with similar 
momenta is due to the bosonic nature of the pions. The phenomenon is often 
referred to as Bose-Einstein Correlation (BEC). The correlation functions carry 
important information about the space-time extent and the dynamics of the particle 
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emitting source. R. Hanbury Brown and R. Q. Twiss used the bosonic correlation 
of particles to measure the angular diameter of distant stars [1]. BEC, or the 
HBT effect has become an important tool in high energy physics. In most heavy 
ion experiments there is enough statistics to measure the two- and three-particle 
correlation function of charged pions. We define the correlation function C(ki, k 2 ) 
as 

N 1 (ki)N 1 (k 2 ) 

Here fci and k 2 are the four-momenta of the emitted particles, Ni(k) is the single 
particle momentum distribution, N 2 (ki, k 2 ) is the two-particle momentum distribu- 
tion. To raise the statistics, one integrates over some components of k\, k 2 . From 
the experimental point of view the adopted definition is 

cm = f^j, (2) 

1 v mix / 

the number of pairs with q relative momentum from the same event divided by 
the number from mixed events. This is measured directly in experiments, and it 
must not to be confused with the ideal Bose-Einstein correlation function without 
any final state effects, Co(q). In our notation q = y— (ki — k 2 ) 2 is the invariant 
relative momentum of the pair, which is equal to 2k i2 in the pairs c.m.s. frame. 
The definition of ki 2 is: k± 2 = \ |ki — k 2 |. 

In the approximate case when multi-particle correlation effects are negligible, 
the BEC has good theoretical description even in the case when final state interac- 
tions (FSI) play important role. The two-particle correlation function of the pions 
can be calculated from the source function as follows [2]: 



C(q)-1 = -J d 3 r(|<f> c ( g ,r)| 2 -l)S(r), (3) 

where $c(<Z,r) is the symmetrized final state two-particle outgoing wave function. 
In case of Coulomb-interaction, this has an analytic form [3,4]. With the imaging 
method (by inverting this integral equation) one can determine the S(r) source 
function. If we assume spherical symmetry, then S(r) depends only on r := |r|. 



2. PHENIX Source Imaging 

Emission source functions have been extracted from correlation functions constructed 
from charged pions produced at mid-rapidity in Au+Au collisions at ^/snn = 
200 GeV [5]. A model- independent imaging technique developed by Brown and 
Danielewicz [6] has been carried out (See Fig. 1). The source parameters extracted 
from low kr PHENIX data give the indication of a long tail, which differs from that 
of a Gaussian distribution. 

The aim of this talk was to present our investigations carried out with THERMI- 
NATOR and HRC Monte Carlo simulatons, that concluded in possible kinematic 
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Fig. 1. (a) Correlation function, C(q), 
(full circles) for ■k + tx + and iv~-ir~pairs 
(200 MeV < k T < 360 MeV) produced 
in the central Au+Au collisions of 0%-20% 
centrality. Open squares represent the re- 
stored correlation from the imaging tech- 
nique. The dotted line represents the result 
of direct correlation fitting. The open trian- 
gles represent the ID angle-averaged correla- 
tion of the 3D correlation function, (b) and 
(c) show the ID source function S(r) and 
4nr 2 S(r), respectively, from imaging (open 
squares), Spheroid fit to the correlation func- 
tion (full circles), and angle- averaging of the 
3D-Gaussian source function (open trian- 
gles). Systematic errors are less than the 
size of the data points. 
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explanations of the long-range components of the radial source function seen in 
PHENIX. 



3. Hadron Rescattering 

It has been shown [7, 8] that calculations based on a hadronic rescattering model 
agree reasonably well with experimental results on single particle spectra, elliptic 
flow, HBT radii in Au+Au collisions at RHIC energies. The model we use here to 
compare PHENIX data with is Tom Humanic's Hadronic Rescattering Calculations. 
The model deals with the 8 most abundant resonances and simulates their decays 
and rescattering. 

Simulated source of charged pions is discussed in three groups, determined by 
time of creation. The core consists of the primordial particles and the decay products 
of the short lived resonances. The secondary pion sources of long lived resonances 
are referred to as the halo, while the decay products of the w(782) meson make up 
a separate class (see Table 1). 

On Fig. 2, note the excellent agreement between PHENIX data and HRC simu- 
lation in the central, low kx case. The agreement is also acceptable in the full range 
for the peripheral low kx case. However, HRC fails to describe the kx dependence 
of the source. 

Further investigations have shown that the w(782) does not account for the 
heavy tail in the simulation. Core-core pairs themselves have a power-law tail in 
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Table 1. Components of charged pion sources in the HRC simulation as grouped 
in this discussion. 



Group 


Lifetime 


Components in HRC 


core 


< 20 fm/c 


primordial tt ± , p, A, K* 


omega 


w 23 fm/c 


w(782) 


halo 


> 25 fm/c 


A, $,77,77/ 



the S l (r 12 ) distribution. In the HRC code, rescattering goes on until it self-quenches. 
The mean free path increases strongly as the system expands, and rescattering in 
such a time dependent mean free path system corresponds to an anomalous diffusion. 
Anomalous diffusion is known to lead to power-law tailed distribution, as stated by 
a generalization of the Central Limit Theorem [9,10]. 

4. Single Freezeout with THERMINATOR 

THERMINATOR— A Thermal Heavy Ion Generator [11], based on the Cracow 
Single Freezeout Model [12] was used for the following simulations. In this model, 
freezeout of the partonic matter occurs on a single space-time surface, governed 
by universal thermodynamic parameters. Particle distribution is determined by 
the thermodynamical equilibrium. Particles propagate freely and the decay of 385 
resonances is fully implemented. In opposition to the HRC simulation, rescattering 
effects are neglected here. 

4-1. PHENIX source simulated 

THERMINATOR simulation provides a flat rapidity distribution. Several checks 
were done to verify that a realistic distribution and pt spectrum of particles is 
generated. 

Centrality can not be directly set in THERMINATOR, but through the pa- 
rameters Pmax and r, where p m ax stands for a radius of a boost invariant axially 
symmetric cylinder, while r denotes the proper time of a simultaneous freezeout 
and hadronization [13]. Although fitted parameters for the spectra of all four RHIC 
experiments including PHENIX were published, PHENIX p T spectra is not well de- 
scribed by THERMINATOR, being primarily developed for, and tuned to STAR. 
Thus we used the parameters corresponding to STAR spectrum fits. This indicates 
that THERMINATOR might be able to show much better agreement with data 
with further tuning. The central data region — 20% has been reproduced using 
simulated events from the three most central regimes in proportion to their span. 
The 40 — 90% regime was substituted with equal number of events simulated in the 
regimes of 50 - 60% and 70 - 80%. 

A separate set of events were simulated for each centrality bin of 5% for cen- 
tral, and 10% for peripheral collisions. Then these sets were concatenated to make 
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Fig. 2. Pion pair source measured in PHENIX and HRC simulation for central 
(0 - 20% ) 200 MeV < k T < 360 MeV collisions (top left); high k T (480 MeV < 
kr < 600 MeV) from the 20% most central Au+Au midy collisions (top right); 
low kT (200 MeV < k T < 360 MeV ) from peripheral collisions (40-90%; (bottom 
center). Solid dots represent data on each figure. Black crosses are simulation, 
while color online symbols stand for different components of the simulated source 
as explained on the legend. The Gauss fit to S(r < 15 fm) is shown as a continuous 
line. 



up the wider regimes of PHENIX data. THERMINATOR was able to reproduce 
the centrality dependence of the data. It is also capable to describe the kx depen- 
dence to a much higher accuracy than HRC does, although there is a significant 
discrepancy at the higher kr region. (See Fig. 3). 

It must be noted that there is a significant excess in the simulated S(r 12) 
at small r values. Our studies have shown that neither the primordial nor the 
resonance particles are responsible for it alone, but the r — ► peak vanishes if all 
the core particles and the short lived (r > 150McV) resonances are excluded from 
the source (see Fig. 4). This can be explained with the freezeout that occurs in a 
single surface in this model, and the lack of rescattering. 

Further simulations indicate that rescattering and resonance decay create very 
different sources for different particles (7r ± ,K ± and p^). 

A study of direction-dependent PHENIX pion-pair source functions using THER- 
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Fig. 3. Pion pair source measured in PHENIX and THERMINATOR simulation 
for central (0 - 20%) 200 MeV < k T < 360 MeV collisions (top left); high k T 
(480 MeV <k T < 600 MeV) from the 20% most central Au+Au midy collisions 
(top right); low kT (200 MeV < kr < 360 MeV) from peripheral collisions (40 — 
90%,) (bottom center). Solid dots represent data on each figure. Black crosses 
are simulation, while color online symbols stand for different components of the 
simulated source as explained on the legend. 



MINATOR Blast Wave and Single Freezeout and models has been carried out by 
P. Chung et al. [14] and concluded that these models, although they do well in the 
"sideways" and "long" directions, fail in the "outward" direction (in the terms of 
the Bertsch-Pratt parametrization [15]). 

4-2. Discussion of results 

THERMINATOR is a simple model that shows a high power in the description of 
particle spectra of the large heavy ion experiments. Although there is no rescatter- 
ing implemented in this model, it reproduces the PHENIX image with the observed 
heavy power-law tail for the lower kx events with a higher accuracy than HRC — 
although it still slightly overpredicts the tail for the high kr regime, and fails for 
r — > because of its inherent limitations. 

The relatively good agreement including the reproduction of the long-range tail 
can be understood by taking into account that the high number of resonances pro- 
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Fig. 4. r -> behaviour of THERMINATOR simulated tt + tt + and tv~tt~ S(r 12 ) 
for central low kr collisions. Plots from top left to bottom right show all pion 
pairs, all pairs from resonance, all primordial pairs, and all pairs coming from 
the short lived resonances V > 150MeV respectively. Solid dots represent data on 
each figure. Black crosses are simulation, while the color online symbols stand for 
different components of the simulated source as explained on the legend. 



vide us with an almost-continuum distribution of lifetimes, increasing virtually to 
infinity. This results in an approximate power-law-like final distribution [16] just as 
anomalous diffusion does [9]. However, THERMINATOR fails to describe correctly 
the intercept A and the HBT radii R S ide, Riong and R ou t [17]. So it is question- 
able that resonance decays alone were sufficient to describe the experimentally seen 
structure of source as Bialas had proposed. 

5. Conclusions 

Source images of PHENIX data have been studied with two fundamentally different, 
simple models. Both were able to reproduce the observed ID images to a certain 
level, but had their limits. The different kinematic mechanisms give similar shapes, 
that are determined by the same underlying mathematical principle (a generaliza- 
tion of the Central Limit Theorem), and so that any of them can be a potential 
explanation of the measured heavy power-law tails. 
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1. HRC simulations describe low kr PHENIX image, but do not have the right 
k t dependence, as for higher fcy values, the simulated S(ri 2 ) 12 still has the 
large tail, but not the PHENIX data. The reason for the heavy tail is rescat- 
tcring, change of mean free path and thus anomalous diffusion. 

2. THERMINATOR simulations provide a reasonable reproduction of the kr and 
centrality dependent behaviour of the ID PHENIX image. The heavy tail in 
those simulations can be explained by the large number of resonances in that 
model. This produces an effect that is similar to the anomalous diffusion. 

3. However, THERMINATOR fails to describe the direction-dependently parametrized 
PHENIX source functions, suggesting that resonance decays alone are not a 
sufficient explanation of the measured structure of the source. 

4. Measurement of S(r) images for w ± ,K ± and p^ 1 could distinguish between 
rescattering and resonance decay effects. 
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